Fourier-Hankel/Bessel space absolute equilibria of 2D gyrokinetics 
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Two global invariants of two dimensional gyrokinetics are shown to be "rugged" (still conserved by the dynamics) 
concerning both Fourier and Hankel/Bessel Galerkin truncations. The truncations are made to keep only a finite range of 
wavenumber k and the Hankel variable b (or z in the Bessel series). The absolute equilibria are used for the discussion 
of the spectral transfers in the configuration-velocity scale space of kinetic magnetized plasma turbulence. Some 
interesting aspects of recent numerical results, which were not well understood, are explained with more satisfaction. 
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[. INTRODUCTION 

It has long been felt comfortable to sit in the chair of ab- 
solute equilibrium thinking about turbulence of neutral and 
conductive fluids'. The basic idea is that the statistical so- 
lutions of the Galerkin truncated system with a subset of the 
Fourier modes may provide information about the turbulence 
physics such as the energy transfer (or cascade) dynamics 
and final states et al. In some cases the absolute equilibria 
spectra can even coincide with cascade ones 2 ' 3 . Such an ap- 
proach is still proving to be able to give new insights to Hall 
magnetohydrodynamic 4 and gyrokinetic 5 turbulence. The lat- 
ter is the first in the kinetic framework with interesting new 
features which deserves some more introductive discussion 
here as the preparation of our work: 



A. the model 

For kinetic systems, such as the gyrokinetics, there are extra 
dimensions of velocities which should be dealt with appropri- 
ately. Zhu and Hammett 5 study the gyrokinetic equation for 
the electrostatic fluctuations in a slab geometry with uniform 
background magnetic field B = B z reads (for a most up- 
dated historical review of the model even appropriate for non- 
plasma physicists, c.f., Rrommes 6 ) 



dg 
dt 



d(y)R 
dR 



dg 
dR 



TTJ^ = 0, 



(1) 



where g(R,v±) defined at the gyrocenter is a component 
of the fluctuating distribution f,f — exp{— q<p}Fo + g + 
(<p)n + h.o.t., deviating from the maxwellian F . Here 
h.o.t. denotes the higher order terms in the gyrokinetic or- 
dering and (-)r the gyroaverage around R: the average of 
any quantity ^ around a ring, of gyroradius p surrounding 
the gyrocenter R, perpendicular (_L) to the magnetic field 

/*(r)<5(r||-R||)<5[|ri-RJ-p(R)]<A- _ 



direction (||), (W) R = 

/*(r)5(r||-R||),5[|r_ L -R L |. 



-p(R)]d 3 r 



27rp(R) 



-R||)5[|ri-Rx|-p(R)]d 3 r 

[We explain that, through- 



out this article, if not particularly specified, J is used, for con- 
venience, to denote actually the definite integral over the en- 
tire domain of the integration variable(s) - space or velocity. 
Using a Fourier representation ^(r) = ^ k exp(— «k • r)^t, 
and considering a straight magnetic field for simplicity here, 

we have (^)r = X]k ex P(~*' 4 ' ^) Jo{k±p) 1 $'k, where Jo 
is a Bessel function.] The normalization by Plunk et al. 7 



brings conveniences for our discussion: The physical (dimen- 
sional) variables having subscript 'p' is normalized as fol- 
lows: t = t p v th /L, x = x p /pth, V = Up/Pth, V = tPpj^i 

v 3 L 

h = hp-^- and Fo = Fo p vf h /nQ. The equilibrium den- 
sity and temperature of the species of interest are no and Tq\ 
the thermal velocity is u t h = y/To/m; the Larmor radius 
is Pth — Wth/^c where the Larmor (cyclotron) frequency is 
il c = qB/m. This equation is closed by using the gyroki- 
netic quasi-neutrality equation to determine the electrostatic 
potential 



2ir / vdv(g) r — a<p — Tip. 



(2) 



For the plasma in a two dimensional (2D) cyclic box, with v\\ 
being integrated out and ± in v± omitted, the whole system in 
wavenumber space are 

d t g(k, v)-zx P J o(pv)<p(p) • qg(q, v) = (3) 

P+q=k 

and 



£(k)=/3(k) / vdvJ {kv)g(k,v) 



(4) 



Here /3(k) 



2tt 



( r , , and r(x) = h{x 2 )e x is an 

exponentially-scaled modified Bessel function. Iq{x) — 
Jo(ix) and r represents the shielding by the species which 
is treated as having a Boltzmann response of some form [see, 
e.g., Plunk et al. 7 ; anisotropic response model, such as those 
respecting the zonal modes, can also be considered as in Zhu 
and Hammett 5 .] 



B. absolute equilibrium in k — v space 

Fourier Galerkin truncation is defined by setting all Fourier 
modes beyond the wave number set IK — {k : k m i n < \k\ < 
k m ax} C the summation over which will be denoted by 
to be zero. Let us start with the 2D case. As Lee 1 did for 
Euler equation for ideal fluid flow, we first observe that the 
dynamics of the "gas" composed by the real and imaginary 
parts of the Fourier modes, denoted by a, satisfies the Liou- 
ville theorem. Actually = 0, where the dot represents the 
time derivative. Then, we proceed to find the canonical distri- 
bution. The only known rugged (still conserved after Fourier 
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Galerkin truncation) invariants for the 2D gyrokinetic system 
are the (relative) entropy relevant quantity 

G(v) = / ^ 2 = £^(M)i 2 

and the mean effective electrostatic potential "energy" 



d 2 r 



E= + 



2/3(k) 



i^(k)r 



(see, e.g., Refs. H and and references therein) where V is 
the volume (area) of the integration domain. Notice that G(v) 
is a function of v, so we have now one plus a continuum of 
conserved quantities, which is the new feature of the problem, 
compared to earlier work for fluids. Here the constant of mo- 
tion S is formed by introducing ccq and a(v), the "(inverse) 
temperature parameters" as the Lagrangian multipliers when 
maximizing the Gibbs ensemble entropy: 



a(v) |ff (k, v)\ 2 dv + a 27r/3(k) 



vdv Jo(kv)g(k, v) / vdvJo(kv)g* (k, v) 



(5) 



with "*" denoting the complex conjugate. The canoni- 
cal statistics of the Fourier Galerkin truncated system is 
then described by the distribution Z -1 exp{— S}, Z = 
jDaexp{—S}. 

The above is the calculation done by Zhu 8 which is the cor- 
responding continuum limit (concerning velocity) of the re- 
sults with discretized velocities by Zhu and Hammer^. We 
should remark that the electrostatic gyrokinetic equations con- 
stitute an integro-differential system. When the system are 
discretized, we should denote the corresponding solution with 
g to distinguish it from the exact solution, g, of the original 
equations. In general, g(v) ^ g(v). So, to be rigorous, for the 
discretized case in Ref. |we should have replaced g(k, Vi) 
with g(k,Vi), as is the case of distinguishing the solution of 
the discretized equation (say, the difference equation) of a dif- 
ferential equation and the solution of the latter (people how- 
ever often neglect such difference for it is not necessary for 
usual discussions.) 



C. the objective of this paper 

Although the calculation in the last subsection respects all 
possible rugged invariants of the original system, especially 
the results with discretized velocity correspond to the current 
gyrokinetic continuum codes, it is not convenient to discuss 
much the physics in velocity scale space (some finite dis- 
cretization scale however still can be discussed as shown in 
Zhu 8 and has partly been verified numerically 9 .) Especially, 
one could ask what would be the case if some other function 
spaces are used for the expansion of velocity (those in Zhu and 
Hammett 5 are the "top hat" bases as used in the continuum 
code.) This is the question we would address in this article. 



In this paper, we will carry out the calculation of 2D gy- 
rokinetic absolute equilibria in Fourier-Hankel/Bessel spaces 
and study the spectral issues based on the analytical results. 
We will show and compare several absolute equilibria of the 
gyrokinetic system with different treatments. All these abso- 
lute equilibria present some similar features for some aspects 
but also demonstrate different merits. We will then explain 
and comment on some recent numerical results. 



II. ABSOLUTE EQUILIBRIA IN 
FOURIER-HANKEL/BESSEL SPACE 

As said, in the past 5,8 we worked directly with velocity and 
respected all the local quadratic invariants G(v), besides E. 
The results are precise and should be observed numerically. 
However, working directly with v, it is not convenient for 
studying the spectral transfers in velocity scale space, which 
requires introducing transformation of v and using the appro- 
priate global invariants. 



A. Using Fourier series and continuous Hankel 
transform 



It is interesting to note that the integration part in the right 
hand side of the quasi-neutrality condition, Eq. |@), is exactly 
the definition of the Hankel transform which will be denoted 
with a breve," That is, 

0(k) = /?(k) J vdvJ (kv)g(k, v) = /3(k)|(k, k), (6) 

which suggests the application of Hankel transform for v. So, 
further doing Hankel transform for the second argument v (to 
b,) we have 



dg(Kb) 
dt 



= l vdvJ (bv)z x ^J?Jo(jpv)tp(j?) • (7) 

p+q=k 

•q / wdwJo(vw)g(q,w). 



As G(v) is conserved for all v, it is equivalent to 
say that, for any reasonable (test) function T(v), Wt — 
J T{v)G(v)vdv is conserved. In the scale, k — b, space, the 
invariants can be represented as 

B = ^/3(k) / S(b-k)\UKb)\ 2 db (8) 



and 



W T = /E/ 9(Kb)Jo(bv)bdb 



g*(k, b')J (b'v)b'db']T{v)vdv. (9) 



We observe that E is rugged with both Fourier- and Hankel- 
Galerkin truncation (keeping only a subset of k and b), and 
this is the same for W when T = 1. The question is whether 
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we would miss any other rugged invariants with the extra 
Galerkin truncation for b. We believe the answer is "No," 
which can be shown following Kraichnan 10 : Integrating out 
v first in Eq. (O we have 



§{k,b)bdb / "g*{k,b')b'db']J(b,b'), (10) 



k 



where J(b,b') = J Jo{bv)Jo(b'v)T(v)vdv. Now, introduce 
the Galerkin truncation (for however continuous b) formally 
in the same way as for discrete k, that is, forcing g(k, b) = 
for b £ B = {b : b m i n < b < b max }. Then for general 

J(b, b'), as d^gk/^ is generally not zero even for b £ B, 
is not zero. Here, 



Wj 



|(k, b)bdb / |*(k, b')b'db']J(b, b') (11) 



as is the meaning of Fourier-Hankel Galerkin truncation. The 
reason is that for g supported only by B, (= 0) con- 
tains not only but also the integration over b and b' 
of dtg(k, 6)6g*(k, b')b'J(b, b'), which is not zero in general, 
with b' in B but b not. In such case, = only when 

J(b, b') = for one and only one of b and b' is in B: As 
b m in and b max are quite arbitrary, the only nontrivial generic 
function satisfying such a property is the Dirac function, that 
is J(b,b') = S(b — b')/b which corresponds to T(v) = 1. 
The trivial one is J(b, b') = uniformly over the full domain, 
or equivalently T(v) = which is physically not interest- 
ing. It is worthy to point out that when T(v) = S(u — v)/v, 
J(b,b') = J (bu)J (b'u) and that W T = G(u), but for the 
reason just explained G{u) is not a rugged invariant. Actually, 
for given feel, the only solution of J(b, b') = for all b' not 
in B seems to be T(v) = or T(v) = 1. So, we believe E 
and 



w = Y,l ll(M)I 2 M> 

1, J 



are the only rugged invariants with the Galerkin truncation of 
both k and b. 

With the arguments given in the above, we conclude the 
canonical distribution be 



^exp { [ - 7i?E k /3(k) / S(b - k)\UK b)\ 2 db - 

rJ2j\§(k,b)\ 2 bdb]/2}, (12) 



which leads to (\g(k, b)\ 2 ) 



with 5h 



k.b ac- 



7wb+7E/3(k)<5i; ib 

quiring 1 for k = b, and otherwise. The calculation with dis- 
cretization of b, say, that parameterized by uniform lattice size 
Ab, gives (||(k, h)\ 2 ) = r wt ,AHr^(k)s, t , 1 (k falls int0 the 
lattice of b indexed by i^) with T\yAb — > -f\y an d Te — > 7e, 



and, Si k! i — > 5k,b- The spectra densities are then 
J B(k,6) = |/3(k) ( 5(6~fc)(||(k,fe)| 2 ) 

_ 7T (3{k)6{b - k) 
2<y w b + j E P(k)S kib 

W(k,b) = b(\Uk,b)\ 2 ) = 



Jwb + 7B / 3(k)(5fc i 6 



(13) 
(14) 



The mean energy is then 



(E)=Y.. JdbE(k,b) = Y,/0',l'), (15) 



with 



£(k,k) 



p(k) 



2 iwk + 7£/3(k) 



(16) 



for k < k max and k < b max ; £(k, k) is zero for k > b max or 
k > k max and is £(k, b max ) - the right limit - for k = b max 

and k < k max , although the value of J_^ x 5(b — b max )db 
is undefined (one however could try to define it through some 
particular nascent delta function.) 



B. Using Fourier and Bessel series 

As in computer simulations and in finite experiments, the 
velocity is finite, it may also be useful to transform velocity 
variable v to Hankel space with Bessel series by taking v to 
be bounded. Suppose velocity is bounded by V, we then have 

g{k,v) = Y,^ 2 [Ji^T 2 U^z)Mzv/V) 

Z 

= E" 7 «(")$(M), (17) 

z 

with z being the zeros of Jo; and, from the orthogonality rela- 
tionship 



we have 



dg(k, z n ) 
dt 



V 



p+q=k 



m 

and the quasi-neutrality condition 

2 1 

vdvj (kv) E "^2 j^y^( k ' z ) J ° 

= Y,B(k,z)Uk,z) 



(19) 



(-) 



(20) 



With the similar arguments as in the continuous Hankel 
transform case, we can see that now only W is rugged with 
regard to both Fourier and Bessel Galerkin truncation and 



^ = EJu 2 ^ 2j r 2 (*> 
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The ruggedness of E is lost due to the fact that B(k, z) 
is in general distributed (over z). The ruggedness of E 
would be recovered when B(k, z) — f3(k)8 ZyZ(k) and <p(k) — 

/3(k)g(k, Z( fe j) where Zn.\ is some zero of Jo: This happens 
when k — . That is, the velocity bound V is taken to be 
k dependent, which was initiated courageously by Plunk and 
Tatsunc 11 . In this way, 



& = E^flOllfc Hk))\ 2 = E k . 2 f /3(k)ll(k, z)\ 2 5. 



M 2 



Note that the dynamical equation ( TT~9b needs to be changed 
accordingly. 

Before going further, we digress to comment that the in- 
troduction of wavenumber dependent truncation of velocity 
amplitude V(k) brings convenience as well as subtleties. The 
obvious point is that V(k) can not be normalized uniformly, as 
inappropriately done by Plunk and Tatsunc 11 in their Eq. (5), 
and that it will affect the relation between the spectra of E and 
W [the relation given in Plunk and Tatsuno's 1 1 Eq. (6) is then 
not correct - see below] which is critical for the Fjortoft argu- 
ment tried in Ref. 1 1 . Physically, in natural phenomena, labo- 
ratory experiment or numerical simulations, the upper bound 
of velocity is unknown or uncontrolled, so it is not very clear 
what exactly the relevance of the k dependent upper bound of 
velocity besides the mathematical convenience (though it is 
obvious that velocity fluctuations at various spacial scales are 
different.) Nevertheless, we can try to continue the absolute 
equilibrium calculation to make things clearer. 

The absolute equilbrium distribution for such wavenum- 
ber dependent upper bound of velocity amplitude is ~ 
exp{ — (oceE + awW)/2} which gives the spectral density 
of E and W: 



E(k,z) 



^(k)^ )Z(fc) |3(k,z)| 2 ) 

7r5 z>2(fe) /3(k) 



a B7 r/3(k)5 z , Z(fc) +4a w V- 2 (k)J^ z (z) 



(21) 



with znA = kV(k) and 



W(k,z) 4 (2V- 2 (k)J^ 2 (z)\UKz)\ 2 ) (22) 
4 



iraEp(k)J?(z)6 ZiZ , k) V2(k)+4a 



w 



with the similar considerations below Eq. ( fT6] l. By definition 



(E) 
(W) 



Z^ k Z^ z I W(k,z) 



(23) 



C. on the spectral transfer 

Note that for the finite velocity and discrete Bessel series 
case 

W(k,z (k) ) = lv- 2 (k)Jr 2 (z (fc) )£(k)//3(k), (24) 



with E(k) = E(k, z {k) ). And, from Eqs. <H2), O and <0 
for the continuous Hankel transform case, we have 



W(k,k) 



2k 
*/3(k) 



£{k,k). 



(25) 



Note that these relations are valid in general, not only for the 
absolute equilibrium spectra. The Fjortoft arguments concern- 
ing the constraints of spectral transfers can be carried over, 
mutatis mutandis, as pioneered by Plunk and Tatsunc 11 , but 
with flawed analysis and inappropriate statements as com- 
mented by the author 12 : For example, the k dependence of 
the upper bound V(k) could affect the relations and should 
have been taken account [even though they were discussing 
the large k limit where j3(k) was approximated as constant 
and J\ {zn.)) was approximated as and the constraints 
can not tell the directions of the transfers, among others. 

The balance equation for studying the energy spectral trans- 
fer can be derived from Eqs. (|7]i and ( [TBl l. 



d£(k, k) 
dl 



"0 



p+q=k ' 

(l(p 5 p)l(q ; w)l*(k,fc))z x p.q, (26) 



where 



A(k,p,w) = 1/[2tt J vdvJo(kv)Jo(pv)Jo(wv)] 

is the area of the triangle formed with legs of lengths k, p and 
w: the area is taken to be infinite if the three legs are not appro- 
priate to close a triangle 13 . The right hand side of Eq. d26l i is 
called the transfer rate function T(k) satisfying T(k) = 0. 
One should calculate the flux from the transfer rate function 
to determine the energy flow in the k — b space and in gen- 
eral it is a very difficult mission (but Kraichnani 4 - was able to 
estimate the signs of fluxes for some particular cases of 2D 
Navier-Stokes turbulence.) Similarly is the case for W trans- 
fer analysis, and also for the velocity bounded case. In this 
paper, we don't attempt to estimate the value of sign of the 
transfer rate or flux. Writing the balance equation down, we 
just mean to remark that there are details in the transfer which 
are not all well described by "macroscopic" arguments, such 
as the Fjortoft constraints or the tendency of relaxation to ab- 
solute equilibria. For example, Eq. d26l l shows that the energy 
transfer is accomplished by the interactions of the triangles re- 
lating both the configuration- and velocity-scale spaces with 
diagonal (p — p) modes, which should be respected by any 
physical conjectures. 

The absolute equilibrium spectra [Eq. (IT~6b and similarly 
Eq. (|2T1 >1 is useful for predicting the transfer directions and 
large-scale structure formulation. This is very similar to the 
case of 2D Naver-Stokes turbulence, except that we now have 
the extra dimension of velocity fluctuation scales. Just as 
Kraichnan discussed for 2D fluid turbulence 14 , the negative 
temperature [je in Eq. ( [TBI l and in Eq. d2"Tl il state 
shows the energy will peak at the lowest modes [note that 
/3(k) is a decrease function of k]. The relevance of equilib- 
rium statistical mechanics with negative temperature to large 
scales of 2D turbulence has various supports starting from 
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Onsager and recently from the Miller-Robert theory 16 , and 
also from the perturbation from |D absolute equilibrium by 
L'vov, Pomyalov and Procaccia 2 . As Rraichnan discussed 14 , 
the negative temperature may due to the large ratio of E over 
W: In (plasma) turbulence, the external forcing could inject 
large amount of energy into the system and/or the collisional 
dissipation of W stronger than that of E could lead to such 
large E/W. 

The possibility of anisotropic /3(k), as discussed in Zhu and 
Hammett 5 with the consideration of zonal flows, will not only 
change the final absolute equilibrium spectra but also the re- 
lations between them and that the transfer dynamics in scale 
space. Note that the /3(k), and then its anisotropy, enters in the 
absolute equilibrium spectra, even for W, only when k = p. 



III. DISCUSSION 

The main resulted formulae for the discussion of spectral 
transfers are the absolute equilibrium spectral densities Eqs. 
( 1131141 ) and ( 121122b . The former ones are even simpler and 
cleaner than the latter which are "contaminated", with V^k) 
and J\(z), by the effects of scale dependent truncation of ve- 
locity amplitude. In the large z limit and with V almost uni- 
form (Plunk and Tatsunc 11 actually took such limits without 
justification, especially for the uniform V) the two cases are 
similar, but they could also be drastically different beyond 
these limits. Usually, the theory by taking the upper bound 
to be infinity may be sufficient and convenient, as in general 
the realistic largest velocity is large enough. 

We have seen that different treatments of the system may 
lead to different versions of absolute equilibria. None of them 
is the exact turbulence solution, however close it could be. 
Which one is physical or physically relevant, or how the phys- 
ical relevance could be made is the critical issue. Concern- 
ing the "as simple as possible, but not simpler" quote (gener- 
ally attributed to Einstein,) a remark about the theory follows. 
Though all absolute equilibria may be precise regarding the 
respective expansion and truncation, for studying the spectral 
transfers in both configuration and velocity scale space the 
one with the continuous Hankel transform of velocity seems 
to have been "as simple as possible" and the one with Fourier- 
Bessel series using k dependent upper bound of velocity is 
"simple", while the one with Fourier-Bessel series using a 
uniform upper bound of velocity is "simpler"; the one with- 
out further transformation of velocity in Zhu and Hammett 5 
however is not simple enough for this purpose. Of course, 
different versions of absolute equilibria have their own mer- 
its and comparisons among them may bring more insights. 
For example, the ones working directly with velocity, with or 
without discretization of velocity 8 , can help to quantify the 
numerical discretization (the current continuum codes use ex- 
actly the same discretizations as in the analysis) effects and its 
possible physical consequences; and, the ones with or without 
upper bound of velocities as presented here can help to quan- 
tify the effects of finiteness of velocities. Actually, turbulence 
is so complicated, so we need various exact (statistical) solu- 
tions of relevant systems (derived from the original system) 



to gain more insights into it. For example, it may be a good 
strategy to start from the theory respecting all invariants of the 
collisionless (inviscid) system 16 for a better knowledge of the 
possible physical relevance of the dropped ones. 

Different versions of absolute equilibria have their common 
features, such as the negative temperature states with energy 
condensation at lowest modes (which should not be very sur- 
prising as the 2D gyrokinetics plasma should reach the hy- 
drodynamic limit, where such features have been well docu- 
mented, in the cold ion limit,) but also differences. How the 
differences of the spectra behavior would be physically rel- 
evant to predict particular features of real turbulence is sub- 
tle and to be examined. For example, the author 8 discussed 
the effect of finite discretization might lead to also conden- 
sation of G(v) at lowest modes, which seems to have been 
confirmed 9 by numerical simulation of the absolute equilib- 
rium ensemble but however still needs to be further examined 
in real ("real" only numerically with collision operators and 
possibly also forcing) turbulence. It is unclear what the phys- 
ical relevance (if any) is when the upper bound of velocity 
truncation is brutally taken to be uniform so that the rugged- 
ness of E is lost (the absolute equilibria simply correspond to 
the case with ole = 0.) 

In three dimension (3D) dynamics, parallel electric field ac- 
celeration, parallel advection et al., besides the slow gyroav- 
eraged E x B drift, take effect. The physical mechanisms of 
spectral transfer are a lot more. In the standard gyrokinetic 
ordering, other parallel nonlinearity are of higher order and 
the parallel dynamics adds only linear terms. While the un- 
changed nonlinearity preserves E and Wt separately, the par- 
allel motions enter to combine them into a single invariant 5 . 
Extra techniques are necessary to extract the velocity scales 
of the mixed linear and nonlinear phase mixing for investi- 
gating into the spectral transfer, though insights of transfers 
in configuration space scales are relatively easier—. It is noted 
earlier 8 that the 3D absolute equilibrium seems to be more 
robust concerning the finite velocity discretization (quantiza- 
tion), we wonder whether this could be the case when the 
quantization is introduced with finite amplitude of the upper 
bound V. 

Finally, the Fjortoft argument as generalized by Plunk and 
Tatsunc 1 1 constraints the possible transfer directions of one of 
the variables when the other one's direction is already known, 
but it does not tell where the system should start and continue 
going. The absolute equilibria, containing already the infor- 
mation of the constraints from the Fjortoft argument, tells the 
system the arrow of time. So, the absolute equilibria analysis 
may be ready to explain many aspects of the spectral transfer 
results of the current numerical simulations, such as those in 
Ref. QjJ. For example, the propagation directions of the modes 
are explained with the tendency of relaxation to the absolute 
equilibria, especially before the regime when and where col- 
lision and/or forcing (if exists) operator takes effect. And, 
especially, our Eq. ( TT~4-b shows that W(k,b) is symmetric 
about the line k — b and then X)|k|=fe b) should have 
accumulate more on the k axis side where more k modes sit, 
which should persist in the nonequilibrium spectrum without 
particular non-symmetric stirring: This simple result seems 
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to explain the non-symmetric behavior in Fig. 6 of Tatsuno 
et al»ii where larger values are on the k axis side (the same 
phenomena has been found by Watanabe. 9 ) We expect that, 

if Eik|=fc ^(M)/^)' with N ( k ) bein 8 the number of 
modes on the k shell, is plotted, the figure should be sym- 
metric about the line k = b: This is not easily seen from the 
transfer rate function in the corresponding balance equation, 
but we can imagine that, as for any symmetric truncations of k 
and b the corresponding absolute equilibria is symmetric, the 
transfer should also be symmetric. 
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